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УДК 681.5 

ОПРЕДЕЛЕНИЕ ГРАНИЦ ОБЛАСТЕЙ ЛОКАЛИЗАЦИИ НУЛЕЙ И ПОЛЮСОВ СИСТЕМЫ 
С ИНТЕРВАЛЬНЫМИ ПАРАМЕТРАМИ 

О. С. Вадутов, С.А. Гайворонский 


Предлагается многопараметрический интервальный аналог метода корневого годографа для исследования динамических 
свойств систем с несколькими интервальными параметрами, линейно входящими в коэффициенты передаточных функций. На 
основе установленных свойств отображения ребер параметрического многогранника на корневую плоскость разработан алго- 
ритм граничной реберной маршрутизации, позволяющий определять образы границ областей локализации нулей и полюсов сис- 
темы. Приводится числовой пример. 


Введение 


1. Постановка задачи 


Одной из основных задач, решаемых при про- 
ектировании линейной системы автоматического 
управления, является получение требуемого каче- 
ства переходных процессов, которое определяется 
полюсами и нулями передаточной функции замк- 
нутой системы. Их желаемое расположение на 
комплексной плоскости обеспечивается при син- 
тезе системы выбором соответствующих настроек 
линейного регулятора. Однако в реальных систе- 
мах физические параметры элементов, как прави- 
ло, точно неизвестны или могут изменяться по за- 
ранее неизвестным законам в определенных диапа- 
зонах. В результате отличия фактических значений 
параметров от номинальных (принятых при синте- 
зе регулятора) будут отличаться от желаемых полю- 
сы и нули системы. Поэтому возникает необходи- 
мость оценки влияния интервальной неопределен- 
ности параметров на динамические свойства сис- 
темы. 

Возможным путем решения данной задачи мо- 
жет быть построение границ областей локализации 
полюсов и нулей системы. Эти области являются 
отображениями на комплексную плоскость корней 
параметрических многогранников соответственно 
характеристического полинома и полинома числи- 
теля передаточной функции замкнутой системы. 
Важную роль при отображении играет характер 
вхождения интервальных параметров в указанные 
полиномы [1]. Здесь мы ограничимся рассмотрени- 
ем систем, интервальные физические параметры 
которых входят в коэффициенты полиномов ли- 
нейно. 

В этом случае для построения границ областей 
локализации корней полинома может быть ис- 
пользована реберная теорема [ 2 ], следуя которой 
необходимо отобразить на корневую плоскость все 
ребра параметрического многогранника. Однако 
рассмотренные примеры применения данной тео- 
ремы показывают, что отображение всех ребер яв- 
ляется сверхдостаточным (образы только некото- 
рых ребер определяют искомые границы). Поэтому 
для уменьшения вычислительных трудностей 
представляет интерес исследование возможности 
нахождения у многогранника интервальных пара- 
метров этих существенных ребер. 


Пусть передаточная функция замкнутой систе- 
мы по задающему или возмущающему воздействи- 
ям имеет вид 

- РШ 



где Р($) и 0(.ч) - полиномы, в коэффициенты кото- 
рых линейно входят интервальные параметры сис- 
темы. Решение поставленной задачи рассмотрим 
на примере построения областей локализации по- 
люсов системы, считая, что для нулей задача реша- 
ется аналогичным образом. Запишем характерис- 
тическое уравнение в следующем виде 

т 

од = 27г4(*)+ВД=о> (1) 

1=1 

где 7 ) - ишервальные параметры, 7 ) < 7 ) < 7), 
7 ] = 7 } тіп , Т[ = 7 /тах ; А ( (з) и 7?Н) - полиномы по сте- 
пеням 5 . 

Так как т параметров заданы своими граничны- 
ми значениями, то параметрический многогран- 
ник, внутри которого Т і могут изменяться произ- 
вольным образом, представляет собой прямоуголь- 
ный гиперпараллелепипед 



содержащий 2"' вершин. Координаты любой точки 
Р т относительно вершины Ѵ д , < 7 = 1,2" определяют- 
ся выражениями 


Т,=Т? +АТ П / = 1, /и; 


( 2 ) 


(т; -т/)< дт; <(т; -т?), ( з) 

где Д Ті - приращение /-го интервального парамет- 
ра, 77 - его значение в вершине Ѵ д . 

Ставится задача: на основе анализа отображе- 
ния ребер на плоскость корней опреде-лить прооб- 
разы границ областей локализации корней уравне- 
ния (1) при любых значениях Т„ і = 1 ,т из области 
Р Т . 


2. Свойства отображения ребер 

Для отображения <р: Р т — > 5, где 5 - множество 
корней ( 1 ), необходимо иметь соотношение, свя- 
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зывающее координаты точек Р т с п корнями ха- 
рактеристического уравнения. Такое соотношение 
может быть получено в результате подстановки в 
(1) выражения (2) 

^ / (^) + АГ 1 -4(5) + 

+АТ 2 ■ А 2 (з) + ... + АТ т ■ Д„(А) = 0, (4) 

т 

где 7У У ( 5 ) = Тр • А і ( 5 ) + В(з) - вершинный 

І = 1 

характеристический полином. 

Введем в рассмотрение ребра Р ^ которые обоз- 
начим 7?/, где / - индекс АТ { , / = 1,/и, <7 - индекс Ѵ д , 
из которой по ребру изменяется Т і . На основании 
(4) запишем уравнение отображения У/ на компле- 
ксную плоскость корней 

В«(х) + АТ і -А і (,) = 0. ( 5 ) 

Пусть (5) является характеристическим уравне- 
нием системы с единичной обратной связью. Тогда 
передаточная функция системы в разомкнутом 
состоянии может иметь вид 

Щ Ч (А Т і , 8)= АТ, ’ Аі(<8 \ (6) 

' ' /9 ? (У) 

Анализируя (5) и (6) с позиции теории корнево- 
го годографа [3] , заметим, что при изменении Д Ті в 
интервале (3) корни (5), стремясь от полюсов 
функции (6) к ее нулям, образуют однопараметри- 
ческий интервальный корневой годограф. Его вет- 
ви назовем реберными ветвями (обозначим 7?У/), а 
их начала и концы - корневыми узлами Тогда 
для отображения <р: Р г — > 5 будут справедливы вы- 
ражения: (р{ Рр) = К8р, <р(ЕД = Ь д . 

По определению многогранника Р т любая его 
грань является прямоугольником на плоскости из- 
менения двух интервальных параметров Т, и 7} из 
одной вершины Ѵ д . В соответствии с этим обозна- 
чим грань С- 1 , а ее образ - 08 На основании (1) 
запишем уравнение отображения плоскости О 4 . 
при снятии ограничений (2), (3): 

Т г Л,(*) + ТуЛ,Ь) + 

+Хг;-4(ц+ад=о. (?) 

к 

Пусть один из корней уравнения (7) имеет вид 
5 Г = а + і(3, г е \,п. Подставляя в (7) координаты 
и выделяя вещественную и мнимую части, получа- 
ем систему двух линейных уравнений с двумя пере- 
менными Ті и Ту 


Т { ■ Ке А 1 (а, /3) + Ті ■ Ке А } (а, /3) + 


+ Ке 


^Т«-А к (а,Р) + В(а,Р) 


= 0 ; 


Т { ■ Іт А 1 (а, /3) + Ті ■ Іт А у (а, /3) + 


( 8 ) 


+ Тт 


^Тр-А к (а,р) + В(а,р) 


0. 


Для уравнений системы (8) характерны следую- 
щие два случая. 

1. Уравнения независимы, и система имеет 
единственное решение 7) = 7” 7} = Т-. Следова- 
тельно, (р'О г ) = Р\ Р = (Т’,Т-), причем точка Р' 
принадлежит плоскости О... 

2. Уравнения отличаются постоянным множите- 
лем. Из этого следует, что в плоскости О с ‘. существу- 
ет прямая г, описываемая любым из уравнений сис- 
темы (8), причем у'Ог) = I- 

Учитывая эти особенности отображения корня 
5,. на плоскость интервальных параметров, рассмот- 
рим границы области 8 Г его локализации при усло- 
вии 

ч>-'(8 г ) = в;. (9) 

Пусть г/г '(.ѵ г ) = Р\ причем Р* е Р/. Так как коор- 
динаты Р* являются единственным решением (8), 
то РУ/ - единственная ветвь, проходящая через я г 
Следовательно, при условии (9) границами 5,. явля- 
ются непересекающиеся образ ы ре бер С 

Пусть (р'Ь,) = I и тс 1 ! = ( точки Р / и Р 7 
принадлежат ребрам Су 7 ). Следовательно, <р(Р 1 Р 2 ) 
= 5,. и через (назовем его особым корневым узлом 
и обозначим II’) проходит множество ветвей кор- 
невых годографов по интервальным параметрам, 
образующим пересекаемые прямой 1 ребра. Угол 
входа каждой такой ветви в ІГ равен углу выхода, и 
поэтому все пересекающиеся в ІГ ветви лежат меж- 
ду двумя пересекающимися там же реберными вет- 
вями. Следовательно, в этом случае при условии (9) 
в состав границы У,, входят пересекающиеся образы 
двух ребер С|. 

Назовем ребра, отображающиеся на границы 
корневых областей, граничными ребрами (обозна- 
чим ОКр), их образы - граничными реберными вет- 
вями (0К8р), а связываемые ими корневые узлы - 
граничными корневыми узлами (СИД. 

Покажем, что у Р т существует единственный 
набор СРр для всех областей локализации компле- 
ксных корней. Пусть (р(Ѵ д ) = 5], 8] = ОС] + )Р Ъ и 
ср(Ѵ д ) = 8 2 , х 2 = ос 2 +іР 2 , причем является гранич- 
ным для У 1 при (р\ У ,) = Р г , а з 2 - внутренним для 8 2 
при ср \ 8 2 ) = Р т . Допустим, что один из интерваль- 
ных параметров 7)- при изменении из Ѵ д выходит за 
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свою границу, определяемую условиями (2) и (3). В 
результате может оказаться, что а^е $ 2 - 

В этом случае (р'Лр € Р т , а (р'($ 2 ) е что проти- 
воречит существованию у и а 2 одного прообраза. 
Поэтому, если хоть одно отображение Ѵ д является 
СЛ д , то такими же являются и все другие ее отобра- 
жения. Из этого можно сделать вывод, что при 
отображении ср: Р т => 8 у Р т существует единст- 
венный набор СКр составляющий граничный ре- 
берный маршрут. 


3. Условие наличия особых корневых узлов 

На основе проведенного анализа можно заклю- 
чить, что реберные ветви могут входить в состав 
границы как полностью, так и частично, причем 
последнее является результатом их пересечения в 
V. Данный случай необходимо учитывать при гра- 
ничной маршрутизации Р г , в связи с чем представ- 
ляется целесообразным заранее знать о возмож- 
ности наличия II в области 8 Г 

Очевидно, что необходимым условием Л &8 Г 
является присутствие прямых / в плоскостях Р т . За- 
метим, что если в двухмерном пространстве какой- 
либо грани существует прямая /, то и в плоскости 
любого сечения Р т , параллельного рассматривае- 
мой грани, также будет соответствующая прямая I. 
Поэтому необходимым условием II е 8 Г является 
наличие прямой I хотя бы в одной из плоскостей 
Р т , имеющих одну общую вершину. 

Как было указано ранее, условием существова- 
ния прямой ( в пространстве параметров Т і и 7^ и ее 
отображения в 7Г(а;у'Д) является линейная зависи- 
мость уравнений (8). Эта зависимость имеет место, 
если существуют такие а и Д что выполняются ра- 
венства 

КеД(а,/ 3) А ^ (а, Р) _ 

Іт А.(а, Р) Іт А/ а, /]) 


Ке 


Іт 


X Т*-А к (а,Р) + В(а,р ) 


X Т<-А к (а,Р) + В(а,р ) 


( 10 ) 


Для этой проверки необходимо решить полу- 
ченную из (10) систему двух нелинейных уравне- 
ний 

[ Ке А і (а, Р ) • Іт А ] (а, (і) - 
-Ке А^а,Р) -Іт А і (а, Р) =0; 

К еА.(а,Р)х 


хіт 


X Т<-А к (а,Р) + В(а,р ) 


(П) 


-Іт А.(а,Р)х 


хКе 


X Т<-А к (а,Р) + В(а,р ) 


= 0 . 


Левые части каждого из уравнений системы (11) 
являются полиномами нечетной степени Р без сво- 
бодного члена. Поэтому, если (11) не имеет реше- 
ния с Р Ф 0, то в отображениях рассматриваемой 
грани и параллельных ей сечений Р т нет Л. Отсу- 
тствие решений (11) для всех сочетаний интерваль- 
ных параметров означает, что границы 8 Г состоят 
из непересекающихся реберных ветвей. 

4. Основные фазовые соотношения 

Определим условия принадлежности Л д грани- 
це 8 Г Для этого введем в рассмотрение угол выхода 
К? из комплексного II д , являющегося полюсом 
функции (6). Обозначим этот угол Ѳ/. Так как на 
функцию (6) распространяются все известные 
свойства корневых годографов, то Ѳ/ можно найти 
из уравнения фаз [3], записанного для Л д . Если 
7>(Д имеет степень п, а А г {а) степень г, тогда Ѳ/ при 
увеличении Т і находится по формуле 

ѳ< = т°-±ѳ,+2®,, 

к = 1 1=1 

а при уменьшении 7} 

к = 1 1=1 

где Ѳ ^ и Ѳ / - углы между вещественной осью и век- 
торами, направленными из Л д соответственно к 
к - му полюсу и к /-му нулю функции (6). 

Проанализируем возможные направления дви- 
жения $ г из Л д , рассматривая области отображения 
граней Р Т с общей вершиной Ѵ д . В каждой из них 
направление^ движения определяется вектором 
Еі ,■ = Е і + Е , где Е і и Е - реберные векторы, за- 
даваемые Ѳ/ и Ѳ/ [4]. Согласно правилам вектор- 
ной алгебры, _Е, ѵ выходит из Л д и лежит внутри об- 
разованного Е. и Е. угла у, причем 0°< х Ё,у<180°. 

Так как из СЛ д выходят две граничные реберные 
ветви, то соответствующие им векторы образуют 
граничный угол (ЕР д , также лежащий в диапазоне 
[0. . . 1 80°). В этом случае все остальные углы Ч*,-,, и, 
следовательно, определяющие их реберные векто- 
ры должны принадлежать (ЕР Выразим данное 
условие через углы выхода реберных ветвей, отсчи- 
тываемые от положительной вещественной полуо- 
си. Оно соответствует выполнению С т 1 неравенств 

|ѳ? - Ѳ?| < 180° . (12) 

Таким образом, на основании проверки (12) 
можно установить принадлежность вершины Р Т 
границе области 8 Г 

Проанализируем возможные соотношения уг- 
лов выхода двух реберных ветвей ( К8, и Р8р из со- 
седних граничных корневых узлов 8 п связанных 
одной граничной реберной ветвью СК8 к . Очевид- 
но, что прообразы Я8, и Я8 І являются ребрами од- 
ной грани Р т . 

Пусть Л" г 8 Г При изменении Т к по ребру СК к 
корень 5,. движется по СК8 к . При этом могут изме- 


66 



Естественные науки 


няться углы выхода из у. ветвей КО/ и КО/. Однако, 
так как ІГ і О г , то КО/ и КО/ не пересекаются. Поэ- 
тому при переходе а г от одного ОІІ д к другому будет 
сохраняться последовательность величин Ѳ/ и Ѳ/. 
Обобщая этот случай на большее число рассматри- 
ваемых реберных ветвей, можно заключить, что в 
соседних граничных корневых узлах значения их 
углов выхода будут располагаться в одинаковой 
последовательности . 

Пусть ІГ е О п тогда КО/ и КО/ могут пересекать- 
ся в ІГ. Если при изменении ^возникнет ситуация 
и г = ІГ , и далее (7* выйдет из О п то на оставшемся 
участке СКО к изменится соотношение Ѳ/ и Ѳ/. В 
этом случае в соседних граничных корневых узлах 
не будет сохраняться последовательность углов вы- 
хода реберных ветвей. 

5. Определение граничного реберного маршрута 

На основе установленных фазовых соотноше- 
ний для граничных реберных ветвей и анализа воз- 
можности их пересечения в особых корневых узлах 
разработаем алгоритм граничной реберной марш- 
рутизации Р т . Пусть системы (11) не имеют реше- 
ний ни для одной из граней Р т с выбранной общей 
вершиной (ІГ ^ 51 г ). Тогда, если известен ОІІ д и вы- 
ходящая из него О ГО Г то очередной граничной ре- 
берной ветвью будет та, модуль разности между уг- 
лом выхода которой и углом выхода известной 
ОКО,' 1 будет наименьшим (наи боль шим). Поэтому, 
зная величины углов Ѳ/, / = 1 ,т для любого гра- 
ничного корневого узла и расположив их в порядке 
возрастания (убывания), можно опре дел ить п осле - 
довательность ветвей КО?, # = 1,2т, /= 1, т, огра- 
ничивающих область локализации комплексного 
корня. 

Этой последовательности соответствует упоря- 
доченный набор ребер Р Т - замкнутый граничный 
реберный маршрут. Логично заключить, что он 
состоит из 2т ребер, связывающих т пар вершин 
Р Т , причем координаты каждой пары имеют проти- 
воположные граничные значения интервальных 
параметров. 

Если хотя бы для одной пары интервальных па- 
раметров система (11) имеет решение, то в области 
0 Г возможно наличие ІГ. В этом случае необходимо 
также найти узел С\] д , затем для него определить 
последовательность Ѳ/ и соответствующую ей пос- 
ледовательность ребер Р ІП . Если при движении по 
полученному таким образом реберному маршруту 
очередное ребро и следующее за ним окажутся об- 
разоваными интервальными параметрами Т ( и 7), 
для которых система (11) имеет решение, то в гра- 
ничный реберный маршрут в этом случае следует 
включить все ребра грани Оу. 

Таким образом, алгоритм определения границ 
областей локализации корней полинома с интер- 
вальными параметрами предусматривает выполне- 
ние следующих этапов: 

1. Приведение характеристического полинома 

системы к виду (1). 


2. Определение координат вершин Р т , соответ- 
ствующих граничным значениям интервальных 
параметров. 

3. Решение систем уравнений ( 1 1 ) в одной из вер- 
шин Р т для установления интервальных пара- 
метров, граничные реберные ветви которых мо- 
гут пересекаться. 

4. Нахождение на основании (12) граничного 
комплексного узла области локализации кор- 
ней. 

5. Граничная реберная маршрутизация Р т по уг- 
лам выхода реберных ветвей из граничного кор- 
невого узла с учетом результатов п. 3. 

6. Отображение на плоскость корней граничного 
реберного маршрута Р т . 

б. Пример 

Для системы автоматической стабилизации си- 
лы натяжения подъемного механизма определим 
(см. рисунок) границы областей локализации кор- 
ней характеристического полинома 
а 2 + а 2 Г + + а 0 = О, 

где я 3 = м//, а 2 = + тхг 2 (\ + Тк І к 2 ), ^ = /с + 

тег 2 + тк \к 2 г 2 (х + Тс), а 0 = тск л к 2 , т - масса гру- 
за, / - длина упругого элемента, /- момент инерции 
электропривода системы, х - удельный коэффици- 
ент демпфирования упругого элемента, с - удель- 
ная жесткость упругого элемента, г - радиус при- 
водного шкива электропривода, к х к 2 - соот- 
ветственно коэффициент усиления электрической 
части привода и коэффициент передачи регулято- 
ра, Т - постоянная времени регулятора. 
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Постоянные параметры объекта управления 
имеют значения: / = 0,5 кгм 2 , % = МО 4 нс, 
с =2- ІО 4 н, г = 0,1 м. Его интервальные параметры 
заданы диапазонами: т е [50,500] кг, / е [50,100] м, 
к 1 е [5,10]. Настройки регулятора к 2 = 1, Т= 0,01 с 
определены для средних (номинальных) значений 
интервальных параметров. 

Приведем характеристический полином систе- 
мы к виду (1): 

I • А 1 ( 5 ) + — • А 2 ( 5 ) + к х ■ А 3 ( 5 ) + А 4 ( 5 ) = 0 , 
т 

где И ; (я) = /У; А 2 (и) = /(с + 

НА = ск 2 г 1+ (ск 2 Тг 2 + хк 2 Н + хТк 2 гѴ\ 

Л 4 («) = (с + 

При трех интервальных параметрах Р 3 содержит 
8 вершин: ^(50,50,5), К 2 (500,50,5), К 3 (500,50,10), 
К 4 (50,50,10), У 5 (50, 100,5), К 6 (500, 100,5), 

Е 7 (500,100,10), К 8 (50,100,10), где первая координа- 
та - т, вторая - /, третья - к х . В результате решения 
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для первой вершины систем (11) установлено, что 
области локализации корней не содержат особых 
корневых узлов. 

Для трех функций вида (6) в первой вершине 
определены постоянные полюсы (3,34+у'3,88, 
-3,34-у3,88, -1,53), а также три набора нулей (три 
нулевых корня МДя); корни Н 2 (у): 0 и -2; корни 
И 3 (я): -2 и -100). Для корневого узла 3,34+/3,88 най- 
дены Ѳ,,, 1 = 34°, Ѳ/ = 5°, Ѳ к / = 86°. При таких углах 
условие (12) выполняется и, следовательно, Ѵ 1 
принадлежит реберному маршруту. Так как, 
Ѳ/<Ѳ,;<Ѳ,Л то имеем последовательность изме- 
нения параметров системы из вершины Ѵ{. I — » т 
— » к х —>/—>• т — » к х —> I. Такая очередность соотве- 
тствует граничному реберному маршруту: Ѵ х — > Ѵ 5 
— > Ѵ 6 — > Ѵ 7 — > Ѵ 3 — > Ѵ 4 — > Ѵ х . Отображение марш- 
рута на верхнюю полуплоскость корней представ- 
лено на рисунке, где построены и все остальные 
несущественные реберные ветви Р 3 . 
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